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Abstract 

It is well known that the ground state of low-dimensional antiferromagnets deviates from 
Neel states due to strong quantum fluctuations. Even in the presence of long-range order, 
those fluctuations produce a substantial reduction of the magnetic moment from its sat- 
uration value. Numerical simulations in anisotropic antiferromagnetic chains suggest that 
quantum fluctuations over Neel order appear in the form of localized reversal of pairs of 
neighboring spins. In this paper, we propose a coherent state representation for the ground 
state to describe the above situation. In the one- dimensional case, our wave function cor- 
responds to a two- mode Glauber state, when the Neel state is used as a reference, while 
the boson fields are associated to coherent flip of spin pairs. The coherence manifests itself 
through the antiferromagnetic long-range order that survives the action of quantum fluctu- 
ations. The present representation is different from the standard zero-point spin wave state, 
and is asymptotically exact in the limit of strong anisotropy. The fermionic version of the 
theory, obtained through the Jordan-Wigner transformation, is also investigated. 


1 Introduction 


The Heisenberg model has been extensively studied for many years as a non trivial many-body 
problem in quantum magnetism. For low dimensional systems, the ground state deviates from Neel 
ordering due to strong quantum fluctuations [1], The determination of this ground state represents 
a fascinating problem, that in one dimension originated a whole branch of Mathematical Physics 
based in the so called “Bethe-Ansatz” technique [2], However, exact solutions are extremely 
intricate, very often not susceptible of a direct physical intuition, and in the case of the Heisenberg 
model, they are restricted to one dimension. 

In this contribution, we would like to present a novel approach based in a localized descriptibn 
of quantum fluctuations. If one takes as a reference the Ising limit, with a ground state of 
Neel type, switching the transverse part of the Heisenberg Hamiltonian may be visualized as 
a disordering process, where pairs of neighboring spins are simultaneously flipped, the ground 
state being a quantum superposition of admixtures contained in the manifold of total S z = 0. 
This effect has been systematically observed in numerical simulations for anisotropic Heisenberg 
chains [3], and was used by Lagos and the author as the heuristic base for the construction of a 
trial solution [1]. To fix ideas, we will concentrate in the case of spin 1 / 2, and most of the examples 
will refer to a one-dimensional system. The theory can be extended to arbitrary dimension [4], 
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and to arbitrary value of the spin [5]. The antiferromagnetic Heisenberg Hamiltonian with axial 
anisotropic exchange can be written as: 

n = J'£ k(>)A(;) + 1 {S + (i)S-(j) + 5-(i)s + 0)}| , (i) 

<ij> L 

where i,j axe site indexes for nearest neighbors, J > 0 is the antiferromagnetic exchange, the S ' s 
axe spin 1/2 operators, and a is the axial anisotropy parameter. Special cases of Hamiltonian (1) 
are: i) a = 0, the Ising case; ii) a — 1, the isotropic Heisenberg model; Hi) a — ► oo, the so called 
‘Xr-modeP. 

It is important to stress that the approach that will be proposed here is not perturbative, in 
spite that the Ising limit is considered as a departure point for its formulation. The method works 
successfully in the axial anisotropic regime, and is asymptotically exact for high anisotropy. Near 
the isotropic point, however, a delocalization transition occurs, and the linear spin wave theory 
becomes a better approximation when compared with exact results or numerical simulations [6]. 
However, the treatment can be extended in a variational way to account for the isotropic case, 
or the Heisenberg-XT regime [7]. In particular, for two dimensions, isotropic exchange, and spin 
S = 1/2, the ground energy deviates less than 0.5% [7] from results obtained by elaborate Monte 
Carlo calculations [8]. 

2 The Wave Function 

The Hamiltonian (1) is said to represent the so called XXZ model, with the axial-anisotropy 
region confined to the interval 0 < a < 1. In the Ising limit (a = 0), the ground state is of Neel 
type. For the infinite chain, there is a broken symmetry, and one of the two possible Neel states 
has to be chosen a s a reference state. They both are connected by time inversion, the ground 
state of the infinite system being a doublet in the anisotropic region. The phase transition, with 
the presence of long-range order and a symmetry broken ground state, requires degeneracy. In 
contrast, in finite chains, the spectrum is not degenerate and the ground state is an eigenstate 
of the time- inversion operator (symmetric or antisymmetric, depending on the number of spins), 
with equal admixtures of both Neel kets. 

For our developments here, we will choose the Neel state |X r ), where the eigenvalues of 
the S z (m) operators for the linear chain are ± (— l) m . With the usual definition of spin ladder 
operators S±{rri), we define boson- like operators by: 

4 = l/f { 7 “^+ E S+(m+l)S-(">)}. (2) 

y iv ^ ** m even J 

4>l = \fjtr\\ aN+ ^ S+(m)S-(m + 1)| , (3) 

where N is the total number of sites in the chain. Operators defined by (2) and (3) flip pairs of 
neighboring spins when applied to the reference Neel state |Af). Two sequences with translational 
symmetry are possible, which we label by even and odd. It is apparent that a similar construction 
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can be realized with the other Neel state | N)' with S z \ (— l) m+1 , interchanging the roles of 
operators (2) and (3). 

In the quasi-Ising limit, the ground state is close to |A/”), and under this assumption we obtain 
the following algebra for the t' s: 

[^e,^] = [ toAl ] = [teAo] = {to At] = 0, (4) 

which are boson-like commutation relations. Within the same approximation, and restricting 
ourselves to the manifold S 2 = 0, the Heisenberg Hamiltonian (1) can be written as a two-mode 
haxmonic oscillator Hamiltonian [1]: 

H = J (M. + fa.) + E,(a), (5) 

where E g (a) is the ground state energy. 

The Neel state satisfies the relations: 

4.W) = \<*\J f W>, hW) = W), (6) 

showing that the Neel state, or quasi-classical state, can be represented as a Glauber state [9], in 
terms of the <f >' s operators. The eigenvalue in (6), that also enters in the definitions (2,3), has 
been chosen so as to cancel the linear terms in Hamiltonian (5). Using a standard notation, we 
define translation operators by: 

D(z) = exp jz</» — z*<f> , (7) 

where <f> may be the even or odd operator. A coherent state is thus obtained as: 

|zi, z 2 ) = exp [zAe ~ z'At} exp \zAo - z* 2 tl} |0), (8) 

with zi,z 2 two arbitrary complex numbers. If we write 

i [n 
z * = z 2 = 2 a V ~2’ 

we get a closed expression for the Neel state a s a minimum uncertainty wave packet of the <£’ s 
operators. Of course, an equivalent representation can be constructed with the Neel ket as the 
vacuum, just by shifting the definitions of the <^’s in a constant, and thus introducing linear terms 
in Hamiltonian (5). 

Since the Neel state is a well defined state, we would like to represent the ground [0 >, in 
terms of fluctuations over the Neel state jA^). This can be accomplished in closed analytic form, 
by inverting expression (8): 

|0) = exp | (<t>i - 4>e) | exp j (<f>l - to) j | AT). (9) 

The structure displayed by (9) is quite interesting. Quantum fluctuations over the quasi-classical 
state |V) are induced by the <^>’s operators. The distribution of fluctuations is Poissonian, the 
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anisotropy parameter a being related to the width of the wave packet. For a sufficiently small, 
the state (9) displays long-range order in spite of quantum fluctuations, but the effective magnetic 
moment is reduced from its saturation value. A vanishing magnetic moment signals at a phase 
transition as a function of the parameter a. 

In spite that the algebra given by (4) is obtained in the Ising limit, the trial state (9) results 
to be extremely accurate in describing the energy and correlation functions in the whole interval 
0 < a < 1. The one dimensional case represents the most stringent test for the wave function, 
since, as we will sketch below, the accuracy of the method improves with the dimension [4]. The 
energy per spin and the staggered magnetic moment corresponding to our trial state (9) can be 
put in closed analytical form in terms of Bessel functions of integer order [5] : 

E g = (0\n\0) = [j 0 2 (2a) + J 1 2 (2a) + 2aJ 1 (2a)] , (10) 

M z = (0|S*(m)|0) = t^J 0 (2a). (11) 

The generalization to higher dimensions is rather straightforward [4]. If one assumes that the 
lattice is bipartite , i.e. not frustrated for Neel ordering, the corresponding boson-like operators are 
defined as: 

(12) 

where R labels sites in a sublattice, {<5} the set of nearest neighbors, and 2 is the coordination 
number. The reference Neel state \J\f >, in this case, assigns up spins to the R sublattice, and 
down spins to R -f <5. The interesting formula is the analogue of (9) for the ground state: 

|0>=exp {"2(I^Tj\/f?W“*)} w - (13) 

Due to the 2 , factor in the exponential of (13), one realizes that the importance of quantum 
fluctuations diminishes with the coordination number 2 , and correspondingly with the dimension- 
ality. However, closed form expression for arbitrary dimension and spin have not been obtained 
yet, even for the simplest lattice structure (like the square lattice) [10]. In spite of the above 
fact, using the Hellmann- Feynman theorem [5], one can show that the error in the energy, as 
calculated with the state (13), is proportional to a 4 . To go towards the isotropic point and to the 
Heisenberg-XT region, one has to generalize the theory in a variational way [7]. 


3 The Jordan-Wigner transformation 

The one-dimensional Heisenberg model for spin 1/2 can be mapped onto a spinless fermion model 
through the so called Jordan-Wigner transformation [11], The two spin states are mapped onto 
fermion states with occupation numbers 0 and 1. Two fermions can not occupy a single site (one 
may think in a very large on-site Coulomb repulsion), and they hop in a lattice with nearest 
neighbors interactions. The transformation is accomplished by: 

S_(n) = exp j— 27 T Y, Cm Cm) C m , (14) 

l m<n J 
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where the C y are fermion operators. We readily obtain the relations: 

5+(m) 5_(m + 1) = C^C m + i, 

S + (m) S-(m) = S x (m) + i = C' m C m = N m , 

that substituted in (1) yield the following fermion Hamiltonian: 

« = t £ (ClCU+1 + * C.) + V £ («. - i) (at« + . - i) , (15) 

where t — Ja/2, V = J. If J > 0, the fermion interaction V is repulsive. For positive a, 
Hamiltonian (15) describes hole- like fermions. An electron-like system is obtained by changing 
a — * —a. We already know that the properties of the Heisenberg Hamiltonian (1) are invariant 
under such a change [12]. Antiferromagnetic ordering implies that the total number of fermions 
is y, and the fermion ket that corresponds to the Neel state is the one with 

W»> = 5(1 +(-!)”]• 

We will call this state \M > to keep continuity with the previous sections. If one translates our 
ground state (9) into fermion language, one gets: 

|0) = exp (16) 

with 

U = 0-1)” (c£ +1 c„ - c^,c ro+1 ) , 

m 

and mean occupation number given by 

Wn|0> = \ (1 + (-l)”Ji.(2a)] , 

where a = 2 t/V. As in the spin case, most physical interesting quantities can be calculated in 
closed analytical form using the trial ket (16). For the fermion-fermion correlation functions, this 
has been done in Ref. [13]. A quantity which is important to investigate possible Fermi liquid 
behavior, is the one-particle momentum distribution [14]. In our case, this can be obtained 
analytically with (16), yielding: 


(Nk) holes = ^ [1 - sin (2a cos A:)] , (17) 

for hole- like fermions, and 

(Nk)part = ^ [1 + sin (2a cos k)] , (18) 

for a particle-like system. The above distributions display a soft variation with the wave number 
k , and therefore no Fermi liquid behavior. The highly correlated limit a — * 0, yields the constant 
value 1/2, with the complete destruction of the Fermi surface of non-interacting particles. 


435 



4 Acknowledgments 

The author acknowledges support from Funda<jao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP, Brazil), and Fundo de Apoio ao Ensino e Pesquisa (FAEP, UNICAMP, Brazil). 

References 

[1] M. Lagos and G. G. Cabrera, Phys. Rev. B 38, 659 (1988). 

[2] H. A. Bethe, Z. Phys. 71 , 205 (1931); R. Orbach, Phys. Rev. 112 , 309 (1958); M. Gaudin, 
La fonction d’ondt de Bethe , (Masson, Paris, 1983). 

[3] D. Medeiros and G. G. Cabrera, Phys. Rev. B 43 , 3703 (1991); ibid. 44 , 848 (1991). 

[4] G. G. Cabrera, M. Lagos, and M. Kiwi, Solid State Commun. 68, 743 (1988). 

[5] D. Gottlieb and M. Lagos, Solid State Commun. 79, 551 (1991). 

[6] G. G. Cabrera, Phys. Rev. B 43, 13476 (1991). 

[7] D. Gottlieb, M. Lagos, and M. Montenegro, Solid State Commun. 81, 729 (1992). 

[8] T. Barnes, D. Konchan, and E. S. Swanson, Phys. Rev. B 39, 4357 (1989). 

[9] R. Glauber, Phys. Rev. 131, 2766 (1963). 

[10] It appears that this problem is related to the formulation of the Wigner-Jordan transformation 
in two dimensions. See for example, Y. R. Wang, Phys. Rev. 46, 151 (1992). 

[11] P. Jordan and E. Wigner, Z. Phys. 47 , 631 (1928); E. Lieb, T. Schultz and D. Mattis, Ann. 
Phys. (N. Y.) 16 , 407 (1961). 

[12] J. des Cloiseaux and M. Gaudin, J. Math. Phys. 7, 1384 (1966). 

[13] M. Lagos, M. Kiwi, E. R. Gagliano, and G. G. Cabrera, Solid State Commun. 67, 225 (1988). 

[14] G. G. Cabrera, to be published. 


436 



